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1.  Introduction 


The  objective  of  this  work  is  to  understand  current  collection  by,  and  the  potential 
and  charge  density  around,  a  high  voltage  probe  in  an  ionospheric  plasma  in  the  absence 
of  scattering.  As  a  canonical  case,  we  have  chosen  a  spherical  probe  10  cm.  in  radius, 
held  at  a  potential  of  11,000  volts,  in  a  plasma  of  density  10  m  ,  and  temperature  0.1 
eV,  with  a  magnetic  field  of  0.3  gauss.  This  plasma  has  a  thermal  current  of  8.46x10^ 
A-m_i  and  a  Debye  length  of  0.743  cm. 

The  best  estimate  for  collected  current  in  this  regime  is  the  Parker-Murphy  limit, 111 
which  we  discuss  in  Section  2.  In  Section  3  we  describe  the  features  of  a  computer  code 
which  we  used  to  simulate  this  problem.  Results  for  the  canonical  case  and  related  cases 
are  described  in  Section  4.  Section  5  discusses  the  results  of  an  orbit-limited  theory  for 
the  charge  density  (and  thus  potential)  about  the  probe.  The  theory  itself  is  derived  in  an 
appendix.  Conclusions  are  given  in  Section  6. 


2.  Parker-Murphy  Limit 

The  Parker-Murphy  limit111  is  derived  by  considering  the  maximum  cylinder  (about 
a  field  line  through  its  center)  from  which  a  probe  can  draw  current,  consistent  with  con¬ 
servation  of  energy  and  canonical  angular  momentum.  The  canonical  angular  momen¬ 
tum  is  given  by 

p0  =  mr2(d0/dt  +  co/2) 

where 

o)  =  eB/m, 

1  Parker.  L  W.  and  B.  L.  Murphy,  "Potential  Buildup  on  an  Electron-Emitting  Ionospheric 
' 1 ’ournal  of  Geophysical  Research  77  >636  H967). 
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and  us  associated  kinetic  energy  by 


(p0  -  mr2(o/2)2/2mr2. 

Neglecting  the  electron’s  kinetic  energy  at  its  maximum  radius,  R,  and  the  r-  and  z- 
components  of  velocity  at  the  probe  radius,  a,  we  find 

e<J>  =  (mR2co/2  -  ma2co/2)2/2ma2 

or 

R2/a2  =  1  +  (8e<))/mto2a2)1A. 

The  maximum  current  to  the  probe  is  then  given  by 

2  7tJlh(R2/a2) 

where 

J,  =  ne  (e0/2jtm)/4. 

th 

For  the  parameters  of  our  canonical  problem  we  have 

to  =  5.275xl06  sec'1 
R  =  1.54  m 
IPM  =  12.6  mA. 


3.  Computational  Method 

We  have  developed  a  2Vz  dimensional  computer  code  (two-dimensional  coordinates 
and  three-dimensional  velocities)  to  simulate  this  problem.  The  code  operates  in  cylindr¬ 
ical  (r,  z,  0)  coordinates  on  a  quadrilateral  grid  consisting  of  90  spherical  shells  extending 
from  the  probe  surface  (0.1  m)  to  a  radius  of  7  meters.  Each  shell  had  25  angular  divi¬ 
sions  ranging  from  the  z-axis  to  the  r-axis.  The  z=0  plane  was  treated  as  a  mirror  plane, 
and  the  r=7m  sphere  was  neld  fixed  at  zero  potential. 

Electrons  were  started  along  the  magnetic  field  (z)  direction  from  just  inside  the  7m 
sphere  with  a  radial  spacing  of  2.3  cm.,  and  with  current  and  velocity  set  to  reproduce  the 
one-sided  plasma  thermal  current  and  one-half  the  ambient  density.  (Beyond  the  collec¬ 
tion  radiu„,  the  remaining  half  of  the  ambient  density  comes  from  electrons  reflected 
from  the  7=0  mirror  plane.)  Energy  and  canonical  angular  momentum  were  explicitly 
conserved  at  each  particle  timestep.  Each  electron  was  followed  until  it  either  (a)  hit  the 
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probe;  (b)  escaped  the  grid;  or  (c)  exceeded  the  limit  on  how  many  times  it  was  allowed 
to  bounce  off  the  mirror  plane.  At  each  timestep  the  particle’s  charge  was  bilinearly 
shared  to  the  nearby  gridpoints.  Ion  density  was  set  to  ne"2^710. 

Iterations  consisted  of  alternately  solving  Poisson’s  equation  and  tracking  particles 
in  the  resulting  potentials  to  update  the  charge  density.  This  process  was  continued  until 
steady  state  solutions  were  achieved.  Potentials  were  constrained  to  be  non-negative,  and 
in  most  runs  charge  densities  were  constrained  to  be  non-positive.  (The  charge-density 
constraint  affects  potentials  along  field  lines  within  the  collection  radius,  which  in  the 
quasi-neutral  limit  satisfy 

e-2<t>/Jie=05(1  +Jt^/0)-Vi< 

It  had  negligible  effect  on  collected  currents  or  on  potentials  beyond  the  collection 
radius.) 

4.  Computational  Results  and  Discussion 

Runs  were  done  for  the  canonical  set  of  parameters  and  for  related  parameters.  A 
complete  list  of  the  runs  with  their  results  is  found  in  Table  I.  Line  5  of  Table  I  is  our 
best  solution  for  the  canonical  case.  The  canonical  case  is  difficult  because  the  Debye 
length  is  short  compared  to  the  mesh  spacing  used,  and  because  the  particle  energy  is  tiny 
compared  with  the  potential  scale  of  the  problem.  Much  insight  can  be  gained  by  calcu¬ 
lating  comparable  cases  with  higher  energy.  Lines  11  and  14  are  our  best  solutions  for 
comparable  cases  in  which  the  plasma  thermal  current  is  the  same  as  the  canonical  case, 
but  the  plasma  temperature  is  raised  to  1  and  10  eV  respectively.  The  other  cases  are 
presented  for  illustrative  purposes. 

The  general  character  of  all  runs  is  that  (a)  particles  originating  within  a  "collection 
radius"  (somewhat  less  than  the  Parker-Murphy  radius)  are  "promptly"  (i.e.,  within  three 
crossings  of  the  z=0  mirror  plane)  collected  by  the  probe;  (b)  most  particles  originating 
between  the  "collection  radius"  and  the  Parker-Murphy  radius  escape  the  sheath  region; 
some  are  collected  after  executing  lengthy  trajectories;  (c)  there  is  a  very  substantial 
region  extending  outward  from  the  Parker-Murphy  radius  whose  particles  are  "trapped", 
i.e.,  escape  the  sheath  only  after  very  lengthy  trajectories  involving  tens  of  crossings  of 
the  mirror  plane;  and  (d)  narticles  traversing  only  the  outermost  portion  of  the  sheath  are 
not  trapped.  How  close  the  "collection  radius"  is  to  the  Parker-Murphy  radius  is 


3 


determined  primarily  by  the  radial  extent  of  the  sheath,  with  'onger-range  potentials  col¬ 
lecting  smaller  currents.  Note  that  Table  I  is  ordered  by  the  radius  at  which  the  potential 
drops  to  100  volts.  The  current  collected  is  very  strongly  inversely  correlated  with  this 
100  volt  radius. 

The  number  of  "bounces"  (i.e.,  passages  through  the  mirror  plane)  allowed  affects 
collected  current  both  directly  (i.e.,  a  particle  that  initially  missed  the  probe  may 
eventually  stumble  into  it),  and  indirectly.  The  indirect  effect  is  that  lc..g  trajectories 
accumulate  a  large  charge  density  which  reduces  the  extent  of  the  sheath.  These  points 
are  well  illustrated  by  the  runs  shown  in  Table  II,  in  which  increased  bounce  number 
leads  to  more  space  charge,  and  thence  to  a  more  compact  potential,  more  collected 
current,  and  less  trapped  current.  Lines  16  and  17  of  Table  I  are  calculations  performed 
using  the  Langmuir-Blodgett  (spherical  space-charge  limited)  potential,  which  happens 
to  be  similar  to  the  potential  of  line  15.  In  going  from  three  bounces  to  50  bounces  (line 
17  to  line  16)  the  collected  current  increases  from  4.6  to  7.6  mA,  still  well  below  the 
Parker-Murphy  limit.  However,  even  after  50  bounces  40  mA  of  current  remains 
trapped.  Taking  into  account  the  effect  of  the  trapped  charge  on  the  potential  (Table  II) 
we  get  a  larger  increase  (to  10  mA)  of  current,  as  the  100-volt  contour  moves  in  from  2.1 
m  to  1 .4  m.  Furthermore,  we  dispose  of  all  but  6  mA  of  the  trapped  current. 

Increasing  the  temperature  reduces  the  collected  current  primarily  through 
decreased  shielding  of  die  potential.  Shielding  is  decreased  because  the  lifetime  of 
trapped  particles  is  reduced,  due  to  increased  probability  of  escape.  Increasing  the 
plasma  temperature  from  0.1  to  1  to  10  eV  (Table  III)  we  find  that  the  100-volt  contour 
moves  outward  from  1.4  to  1.6  to  2.0  meters,  and  the  collected  current  decreases  from  10 
mA  to  9  mA  to  6  mA.  Note  that  the  current  is  reduced  by  the  shape  of  the  potential, 
rather  than  directly  by  the  higher  energy  of  the  incident  particles.  Line  4  of  Table  I 
shows  a  10  eV  plasma  with  ten  times  the  thermal  current  of  the  canonical  case.  Its  poten¬ 
tial  profile  is  comparable  to  that  of  line  6,  and  its  collected  current  scales  directly  with  the 
plasma  thermal  current,  despite  a  hundredfold  increase  in  particle  energy. 


2.  Langmuir,  I.  and  K.  B.  Blodgett,  "Currents  Limited  by  Space  Charge  between  Concentric 
Spheres",  Physical  Review  24, 49-59  (1924). 
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Runs  which  allowed  positive  charge  density  had  positive  potentials  extending  along 
the  field  lines  within  the  collection  radius.  This  had  no  effect  on  the  collected  current. 

Figures  1  and  2  show  the  electrostatic  potentials  and  particle  trajectories  for  the  0.1 
eV  and  10  eV  cases  (lines  5  and  14  of  Table  I).  (Note  that  the  actual  runs  followed  more 
trajectories  for  longer  paths  than  shown  in  the  figures.)  The  positive  potentials  extending 
along  the  field  lines  are  seen  clearly  in  the  10  eV  case.  Note  also  the  line  of  "turning 
points"  for  particles  which  missed  the  sphere.  These  particles  are  actually  orbiting  the 
symmetry  axis,  and  contributing  substantially  to  the  charge  density.  Both  of  these 
features  are  also  present  in  the  0. 1  eV  case,  but  are  difficult  to  see  in  the  figure. 

Figure  3  shows  the  charge  density  for  the  0.1  eV  case  (Table  I,  line  5).  Note  that 
the  charge  density  is  relatively  low  in  the  region  traversed  by  collected  electrons,  and 
higher  in  the  trapping  region.  Figure  4  shows  the  charge  density  along  a  radial  line  nor¬ 
mal  to  the  magnetic  field  for  all  three  temperatures.  For  all  but  the  10  eV  case,  there  is  a 
sharp  peak  in  the  charge  density  at  a  radius  of  about  0.4  m.  This  peak  consists  of  elec¬ 
trons  which  originated  just  outside  the  Parker-Murphy  radius.  For  a  circular  orbit  at  a 
given  radius,  the  centrifugal  force  (outward)  and  the  magnetic  force  (outward)  on  an 
electron  are  determined  by  conservation  of  orbital  angular  momentum.  These  do  indeed 
balance  the  electric  force  (inward)  at  a  radius  of  about  0.4  m.  However,  the  kinetic 
energy  associated  with  this  circular  motion  is  well  under  half  the  total  kinetic  energy  of 
these  electrons;  most  of  their  kinetic  energy  is  associated  with  oscillations  parallel  to  the 
magnetic  field. 

5.  Orbit-Limited  Theory 

The  appendix  presents  a  theory  for  an  upper  bound  to  the  charge  density  for  this 
problem.  By  "orbit-limited"  we  mean  that  any  point  in  phase-space  whose  energy  and 
canonical  momentum  are  consistent  with  a  trajectory  connecting  to  ambient  plasma  con¬ 
tributes  to  the  local  density.  This  is  known  to  give  a  charge  density  which  is  far  too  high, 
because  it  contains  contributions  from  trajectories  which  can  reach  a  field  point  only  in 
the  presence  of  bizarrely  unrealistic  potentials.  However,  we  can  modify  the  orbit- 
limited  theory  as  follows  to  assure  that  the  charge  density  is  built  up  from  only  "reason¬ 
able  '  trajectories. 
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We  first  note  that  trajectories  originating  within  a  "collection  radius",  Rp  are 
immediately  collected  by  the  probe.  Therefore,  we  exclude  this  region  of  angular 
momentum  from  the  orbit-limited  charge  density  integrals,  and  insteaa  treat  such  trajec¬ 
tories  by  accelerated  current  formulas.  We  also  exclude  particles  outside  some  larger 
radius,  R2,  which  may  be  roughly  interpreted  as  the  "range  of  the  potential",  on  the 
grounds  that  such  particles  never  see  any  fields  which  draw  them  into  the  sheath  region. 
Outside  R2  we  use  accelerated  current  formulas  to  determine  electron  density. 

If  we  take  the  "range  of  the  potential"  to  be  the  radius  at  which  the  potential  drops 
to  the  plasma  temperature,  we  can  imagine  solving  the  problem  self-consistently  using 
these  "modified  orbit-limited"  charge  densities.  This  can  be  done  easily  only  for  the  case 
of  the  10  eV  plasma.  We  find  that  the  resultant  potential  is  qualitatively  correct,  but  too 
compact,  i.e.,  the  charge  density  is  still  too  high.  It  follows  that  the  currents  in  this 
potential  will  be  too  high.  Figure  5  shows  the  solution  for  the  10  eV  case  using  the 
modified  orbit-limited  theory. 

Alternatively,  we  can  generate  any  reasonable  set  of  potentials  (either  by  using  a 
particle  tracking  code  or,  for  example,  by  solving  Laplace’s  equation  inside  a  grounded 
ellipse),  set  Rj  to  be  at  or  just  inside  the  Parker-Murphy  radius,  and  determine  R2  by  the 
requirement  that  the  total  charge  in  the  sheath  cancel  the  charge  on  the  sphere.  (This 
method  does  not  require  that  R2  -  Rj  be  resolvable  by  the  computer  code.)  Figure  6 
shows  charge  densities  calculated  for  the  0.1  eV  case  using  the  potentials  of  figure  1,  Rj 
=  1.54  m,  and  R2  =  1.56  m.  Comparing  figure  6  with  figure  3,  we  find  that  the  theory 
gives  a  much  sharper  boundary  between  the  high  and  low  density  regions,  and  that  the 
theory  does  not  predict  the  structure  due  to  particle  dynamics  in  the  high  density  region. 
It  can  be  shown  (see  appendix)  that  the  theory  predicts  a  density 

p/p0  =  (mw2r2/87iT)‘'*  x  [(R^-R^/r2] 

for  all  points  whose  radius  and  potential  are  consistent  with  trajectories  connecting  to  the 
spherical  shell  Rj<  R  <  R2. 

Table  I  includes  results  from  the  "modified  orbit-limited"  theory  for  the  0.1,  1,  and 
10  eV  cases.  Only  in  the  10  eV  case  is  the  answer  truly  self-consistent;  the  charge  densi¬ 
ties  in  the  other  two  cases  are  too  high  to  easily  achieve  self-consistency.  However, 
these  potentials  follow  the  same  theme  of  giving  a  potential  which  is  too  compact  and 
draws  more  current  than  obtained  by  particle-tracking. 
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6.  Conclusions 

We  have  used  a  particle-tracking  code  to  calculate  potentials  and  currents  associ¬ 
ated  with  an  11  kV,  10  cm  probe  in  an  ionospheric  plasma.  We  have  also  developed  a 
"modified  orbit-limited"  theory  which  yields  potentials  about  the  probe  which  are  quali¬ 
tatively  good,  but  quantitatively  are  too  compact. 

For  fixed  probe  voltage  and  magnetic  field,  the  probe  current  has  a  strong  inverse 
correlation  with  the  extent  of  the  potential,  i.e.,  longer-ranged  potentials  draw  less 
current.  The  reason  for  this  is  that  particles  at  a  given  radius  receive  more  longitudinal 
acceleration,  and  less  radial  acceleration,  as  the  potential  becomes  longer-ranged.  Parti¬ 
cles  originating  at  the  Parker-Murphy  radius  can  hit  the  probe  only  with  zero  longitudinal 
velocity,  and  the  probability  of  having  a  small  longitudinal  velocity  when  near  the  probe 
is  less  likely  for  a  longer-ranged  potential. 

The  range  of  the  potential  is  determined  by  particles  which  are  "trapped".  Most  of 
these  originate  in  a  region  just  outside  the  Parker-Murphy  radius,  and  spend  a  large 
amount  of  time  orbiting  close  to  the  probe.  This  is  illustrated  by  figure  4,  which  plots  the 
net  charge  density  for  the  three  cases.  Note  the  very'  high  (several  times  ambient  density) 
peak  in  the  charge  density  at  about  0.4  m.  This  corresponds  approximately  to  the  radius 
of  a  circular  orbit  where  electrostatic  (inward),  centrifugal  (outward),  and  magnetic  (out¬ 
ward)  forces  balance  for  electrons  originating  just  outside  the  Parker-Murphy  radius  and 
moving  in  the  (ExB)  direction. 
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Table  I. 

Potential  and  Current  Collection  Calculations 


# 

V 

IV] 

T 

[eV] 

n 

[m3] 

+Q 

Allowed 

NBounce 

Vo] 

[mA] 

^r«p 

1mA] 

R(100) 

]m] 

R(10) 

[m] 

R(0.1) 

M 

1 

11.000(A) 

0.1 

lxlO11 

n 

23 

11.7 

5 

1.0 

1.45 

1.9 

2 

11.000(B) 

1.0 

3xl010 

n 

23 

12.4 

6 

1.1 

1.7 

2.3 

3 

11.000(B) 

1.0 

3xl010 

y 

- 

- 

- 

1.1 

1.7 

2.4 

4 

11,000 

10. 

lxlO11 

n 

13 

110 

32 

1.3 

1.7 

1.9 

5 

11.000 

0.1 

1x10“ 

y 

50 

iO 

6 

1.4 

1.8 

1.9 

6 

11,000 

0.1 

lxlO11 

n 

50 

11 

9 

1.4 

1.8 

2.0 

7 

11,000 

0.1 

1x10“ 

n 

20 

9.0 

22 

1.5 

1.9 

2.0 

8 

11,000(J) 

10. 

lxlO10 

n 

23 

7.3 

7.2 

1.5 

2.2 

3.2 

9 

11.000(C) 

10. 

lxlO10 

y 

23 

- 

- 

1.7 

2.6 

3.6 

10 

11,000 

1. 

3xl010 

n 

23 

8.0 

8.0 

1.6 

2.1 

2.3 

11 

11,000 

1. 

3xl010 

y 

23 

9.0 

5.0 

1.6 

2.1 

2.4 

12 

11,000 

10. 

lxlO10 

n 

23 

6.5 

3 

1.7 

2.2 

2.9 

13 

11,000 

0.1 

lxlO11 

n 

10 

6.7 

30 

1.8 

2.4 

7.7 

14 

11,000 

10. 

lxlO10 

y 

50 

6.8 

0 

1.95 

2.7 

4.5 

15 

11,000 

0.1 

lxlO11 

n 

3 

4.2 

46 

2.1 

2.8 

3.5 

16 

11.000(D) 

0.1 

1x10“ 

- 

50 

7.6 

40 

2.2 

2.9 

3.5 

17 

11.000(D) 

0.1 

lxlO11 

- 

3 

4.6 

46 

2.2 

2.9 

3.5 

(A)  Modified  orbit-limited  charge  density  ;  (R1,R2)=(  1.40, 1.46). 

(B)  Modified  orbit-limited  charge  density  ;  (R1,R2)=(1.50,2.10). 

(C)  Modified  orbit-limited  charge  density  ;  (R1,R2)=(1.13,2.20). 

(D)  Langmuir-Blodgett  Potential. 


X 


Table  11. 

Effects  of  "bounce  number"  on  the  potential 
and  collected  current  calculated  for  the  probe. 


;  » 

V 

T 

n 

+Q 

NBounce 

^Col 

L, 

l  rap 

R(100) 

R(10) 

R(0.1) 

- 

ivi 

IcV] 

[m  3j 

Allowed 

- 

[mA] 

[mA] 

[m] 

[m] 

[m] 

!  6 

1 1  ,000 

0.1 

lxlO11 

n 

50 

11 

9 

1.4 

1.8 

2.0 

7 

1 1 .000 

0.1 

ixur 

n 

20 

9.0 

oo 

1.5 

1.9 

2.0 

■  13 

1 1,000 

0.1 

1  x  1 01 1 

n 

10 

6.7 

30 

1.8 

2.4 

2.7 

15 

1 !  rxx) 

0.1 

1x10  s  1 

n 

7 

4.2 

46 

2.1 

2.8 

3.5 

Table  III. 

Effects  of  increasing  'emperature  (at  constant 
current  density)  on  the  potential  and  collected  current. 


i  n 

i  ' 

V 

(V] 

T 

feV] 

n 

[m  ’] 

+Q 

Allowed 

NBounce 

*Col 

[mA] 

^Trap 

[mAl 

R(100) 

[m] 

R(10) 

[m] 

R(0.1) 

[m] 

!  5 

11,000 

0.1 

1  x  1 0 1 ! 

y 

50 

10 

6 

1.4 

1.8 

1.9 

1! 

11, (XX) 

1. 

3xl010 

y 

23 

9.0 

5.0 

1.6 

2.1 

2.4 

14 

11,000 

10. 

1x10'° 

y 

50 

6.8 

0 

1.95 

2.7 

4.5 
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Figure  Captions 

Figure  1.  Electrostatic  potential  contours  (solid  lines)  and  particle  trajectories  (dotted 
lines)  for  a  0.1  meter  radius  sphere  at  1 1  kV,  magnetic  field  0.3  gauss,  plasma  den- 
sity  10  m  ',  plasma  temperature  0.1  eV.  Contours  are  logarithmically  spaced, 
with  the  outermost  contour  at  1  volt.  The  figure  has  mirror  symmetry  at  the  R-Axis 
and  rotational  symmetry  about  the  Z-Axis.  Distances  are  noted  in  meters. 

Figure  2.  Electrostatic  potential  contours  (solid  lines)  and  particle  trajectories  (dotted 
lines)  for  a  0.1  meter  radius  sphere  at  11  kV,  magnetic  field  0.3  gauss,  plasma  den- 
sity  10  m,  plasma  temperature  10  eV.  Contours  are  logarithmically  spaced,  with 
the  outermost  contour  at  1  volt.  The  figure  has  mirror  symmetry  at  the  R-Axis  and 
rotational  symmetry  about  the  Z-Axis.  Distances  are  noted  in  meters. 

Figure  3.  Charge  density  for  a  0.1  meter  radius  sphere  at  11  kV.  magnetic  field  0.3 

11  3 

gauss,  plasma  density  10  m  ,  plasma  temperature  0.1  eV.  (Compare  potentials 
and  trajectories,  figure  1.)  The  contours  are  logarithmically  spaced  and  denote  the 
negative  of  the  charge  density  (p/e0);  ambient  is  1807.  Note  the  low  charge  density 
along  the  Z-Axis,  the  sharp  peak  on  the  R-Axis  at  R  =  0.5  meters,  and  the  ridge 
extending  upward  from  the  R-Axis  just  inside  the  sheath. 

Figure  4.  Charge  density  (p/e0)  for  the  0.1  eV  case  (curve  with  two  sharp  peaks,  ambient 
=  1807,  compare  fig.  3),  1.0  eV  case  (curve  with  two  broad  peaks,  ambient  =  571), 
and  the  10  eV  case  (flat  curve,  ambient  =  181). 

Figure  5.  Electrostatic  potentials  and  panicle  trajectories  for  a  0.1  meter  radius  sphere  at 

1 0  3 

11  kV,  magnetic  field  0.3  gauss,  plasma  density  10  m’  ,  plasma  temperature  10 
eV,  calculated  by  the  ''modified  orbit-limited  theory”.  (Compare  with  figure  2.) 

Figure  6.  Charge  densities  calculated  for  the  0.1  eV  case  using  the  potentials  of  figure  1, 
Rj  =  1.54  m,  and  R2  =  1.56  m.  (Compare  with  figure  3.) 
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APPENDIX 


ORBIT-LIMITED  CHARGE  DENSITY 
IN  MAGNETIC  FIELD 
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Orbit-Limited  Charge  Density  in  Magnetic  Field 

With  cylindrical  coordinates,  B-field  along  Z,  the 
Lagrangian  is 

r  m  (•  2  >2  2t2)  \  u  r2  % 

L  =  2  [r  +  z  +  r  <p  j  +  etp( r,z)  +  eB  6 

The  canonical  moments  are  [w  =  U)Q  =  eB/m] 

=  mr^$  +  eBr^/2  =  mr^(§  +  w/2) 

pr  =  mr 

• 

p  =  mz 
z 


The  energy  is 


E  ■ u  +  £3  h 


where  U 


_  i_  f  2 

2m  lPr 


2 

mr  u 


2 


e  <p( r,z) 


dp  dp  =  47TmdU 
r  z 


We  need  to  relate  the  distribution  factor,  f  (r , z , 6 , pr , pz , pg) 
to  the  ordinary  density.  Over  some  region  we  have 


N 


p  (r , z , 6) rdrdzd# 


r  r 

drdzdS  dprdpzdp^  f 


so  that 


r 

fdprdpzdp0  =  r p 

J  * 
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In  undisturbed  plasma,  the  distribution  function  is 


fo(R,zo,e  >Uo,p0)  =  Po(2^ml3] 


-  1/2  -U  /T 


-  (p0  -  mR2w/2)2/2mR2T 

x  exp 

I.  Non-Magnetized  Case 

To  do  the  nonmagnetized  case,  we  set  ui  =  0. 

From  a  point  (r,z,0)  within  the  sheath,  trace  an 
orbit  back  to  a  cylinder  of  radius  R,  where  R  lies  beyond 
the  range  of  the  sheath .  Then 


Lkill 

Rp  (R) 


dU  27rmR  T 


2^31  1/2 


Uo*  ^]/T 


where  the  energy 


E  =  Uq  +  p2/2mR2  =  U  -  e<p  +  p2/2mr2  >  0 


and  E  >  pg/2mR2 


2  2 

Pe  pe 

so  that  U1  =  - ~  ^  +  ei p 

2m R  2mr 


19 


and 


lp.Ll ) 

Rp(R)  ' 


/*> 

dpe 


dL 


x^27rmR  /T  J 


where 


Region  of  integration  (w  =  eB/m  =  0) . 


[2ms2(e^-U)J 1/2 


ALLOWED 


21 


to  t*— • 


"'o  that 


£_ 

P  o 


e<p/l 

e 


1 -erf 


1 

2 


Using  erf  ( z ) 


2 

e  Z  /z  -1 7T 


for  e$/T  >>  1 


[T/Tre^J 


1 

2 


4ei^/7rT 


For  e<p / T  <<  1  ,  we 


get 


£_ 


(e^>/T)  +  0 


( e0/T)  2 


II  Magnetic  Case 


eB 

m 


*  0 


As  before,  we  have 


(2*mT3) 1/2rp(r) 


f00 

/■ 

dp0 

'  oo 

dU  p(R)exp 

1 

;t 

J  -oo 


e^ 


2 

*  £4 

2mr‘ 


(Pg-mr  u/ 2) 
2 


(pfi-mr2w/2) 2 

e^-U-  - - 2 - 

2mr 

2  (Pg-£ER2a//2)2 


2raR2 


Now,  we  need  a  new  criterion  for  R,  the  iadius  at  which  we 
terminate  outwardly  tracked  orbits.  Orbits  with  pg  >  0  do 
not  include  the  r=0  axis.  We  terminate  these  orbits  with 


5=0 


R  =2  pg/mw 


f  2  1 

[ Pg~mR  w / 2 


Orbits  with  pg  <  0  do  include  the  r=0  axis;  for  these 

5  =  -uj 


R  =  -2p0/m u 


f  pg-mR^w/2j 


The  region  of  integration  i s 


e<p  - 


f  2 

!  pfl-mr  uj/2 

N  C7 


(Pfl  >  01 


[pe-mr2w/2] 


-  wp e  (pg  <  0) 


23 


(27rmT3)1/2rp(r)/p(R)  = 


r00 

f00  r 

dU 

dp^exp 

e<p-U 


(Pe~nir2w/2)  2 


2mr‘ 


f  e<p 


dU 


'®§-JS  [1+  (8(e^-U)  /mr2w2)  1/2] 
dpg  exp  [•••] 


2  2/c 
ef-mr  w  /8 


2 

mr  D 


1-^8 (e0-U) /mr2w2] 


a  2  2  ,0  ,o 

e<p  -mr  w  /8  ^ 


dU 


dpfl  exp  [•••] 


2 

mr  uj 


l+(8(e^-U)  /mr\-2) 


'e<p 

dU 

e^-mr2u2/8 


2, 

f  mr  u 


2 

dp0  exp 


2 

mr  uj 


1  + 


(8(e0-U) /mr2w2) 1/2] 

[*••] 

1 

(8(e^U)/mr2u2)2  ] 


First  term  gives 

TPp  ^  =  r  exp(e^/T) 
'o 

Second  term  gives 


/•  e<p-mr  u;  /8 

r 

i  1-, 

r 

2 

dU  e(e^-U)/T 

n  2  r  2  2.  2 

9  ♦  ]  J 

, 

o 
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Third  term  gives 


r  e(e^-U)/T  erf^ 


e<p  - 


2,  2 

mr  to 


Fourth  term  gives 


,  2  2  /Q 
ep-mr  w  /8 

ay  e(^-u)/T 

T 


1_ 

2  2%  2 


-erf 


rmr_a>_i 

l  8T  J 


mr  to 


Fifth  term  gives 


1 

(e0-U)/T  j, T  2  2,2  f 

~ T  erf  l"2f"J  1  + 


2  2 
mr  to 


,  2  2  /Q 
ef-mr  w  /& 


'  ,  2  2.  2  , 

erf  t~2T  )  (J 


mr  to 
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Gathering  terms,  we  have 


=  exp  ( e<p/T ) 
P  o 


■e<p 

f  e(e*-U)/T  erf(efcUj1/2 

o 


1 

+  2 


>  2  2 
retp- mr  a/  /8 

|U  e(e*-U)/T  erf[^^ 


1/2 


1 

2 


■  e<p 

f  e(e0^)/T  erf 


2.  2± 


rmrw}2(,  f2fe0-im21 

l  2T  J  l1  l  22  J  J 


mr  w 


o 


1  1 

1 

+  2 

f  e(e^U>/T  erf 

mr  ui  J 

■4>- 


2  2 
mr_w_ 

8 


The  first  two  terms  given  the  B=0  resalt.  In  the  remaining 
terms,  the  leading  term  of  erf  (1)  cancels. 

Expanding  erf,  we  get 
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Term  3: 


Term  4 : 


_ 1 

2^ 


2  2 

,e<p-m r  u)  /8 

dU 

T 

o 


1/2 


il 


+ 


,e<p 

1  dU 

2-J~7r  J  T 


e 


2  2 


X 


r  mr 


2  2, -1/2 


l  2T 


r) 


M^H'1 

mr  u/ 


2  2  /OT 
-mr  u/  /21 


e^> 


,-2  [  (e<p-U)  (mr2w2)/2TY£ 


,  .  f.2.(elzUll 

1  L  2  2) 

mr  ui 


1/2  T 


=  smal  1 
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Term  5: 


1 

2Jlr 


re<P 


dU 
T  e 


-[^^T  ~  (e<p-\J)  (mr2w2)  /T2]  2j 

[mr2w2/2T)2  -  ((etf>-U)/T]2 


^  mr  u 
’  8 


(steps  omitted) 
1  (  2T 

i —  x  2  2 

41  7T  Imr  ui 


1 

2 


So,  finally, 


£_  = 

^o 


III.  A  More  Reasonable  Orbit-Limit. 

We  can  get  a  more  reasonable  orbit-limited  bound  on 
density  by  noting  that: 

a)  Particles  originating  within  radius  Rj ,  hit  the 
object,  and  contribute  negligibly  to  the  density 

b)  Particles  originating  outside  radius  R2  (the 
range  of  the  potential)  may  be  ignored,  as  they 
never  see  the  potential;  and 

c)  We  also  ignore  orbiting  particles  (pg<0) . 
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Def ine 


=  mR^2  w/2 

P2  =  mR22  u/ 2 

P  -  mr2  w/2  [l  +  [8e<p/mr2w2)1/2] 
m3>x 

P  =  max  (0,mr2w/2  [l  -  (8e^/mr^w2) 
min  \  7  /  Tf 

From  above  we  have 


P  In.  41 

Po 


(27rmr2T3)-1/2 


x 


r  00 

]  dU  exp 


max(O,e0  - 


[  (e^-U-[p£  -  mr2w/2] 
[p0  -  mr2w,/2]2 
2mr2 


Case  A: 


P2>  Pl>  P 


max 


Leads  to 


PA(r,<P) 

Po 


e<p/T 


[  erf 

Pg-  mr2w/2 

-  erf 

f  2  'l 

P^-  mr  w/2 

l(2mr2T)1/2J 

30 


2/2mr2]  /T 


If  the  arguments  of  erf  are  large,  we  may  use  the  asymptotic 
expansion  to  write 


V  mA/2  P2-  Pmay  f  ea-ll/2 
bmAO1/2  '  (i^TJ1/2  *  1  T  J 


Pg-  mr^w/2 

1 

f(2mr2T)1/2j 

l(2mr2T)1/2J 

'  2TT 

P0-  mr2w/2 

^  J 

(P2  Pmax3 

(P2  Pmax  f  1/2 

x  exp 

[*  (2nr2T)1/2 

[(2mr2T)'/2  iT  J  J 

(and  similarly  for  terms) 


Case  B: 


P  >  P„>  P,>  mr2w/2  [l 
max  2  1  L 


8ed>  1  1/2  ' 

2  ,2 

Imr  uj  J 


In  this  case  the  inner  integral  simply  gives  T,  so  we 

PB(r,<P ) 


f  2  2) 
mu  r 

1/2 

[R22-  Rl2' 

l  87rT  J 

2 

l  r  J 

27rmr2T  J 


1/2 


(P2-  Pl) 


Case  C:  mr  Ui/2  [1 


(  8e$  ' 

1/2  ‘ 

2  2 
Imr  w  J 

>  Pr 
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This  is  similar  to  Case  A,  giving 


PC(r^)  1  „  e<p/T  f  f 

mr2w/2  -  Pj 

p  "  2  °  l  f 

ro 

dz-AP'2  J 

r  2  1 

mr  a//2  - 

UzmAO1/2  J 

For  large  argument,  we  have 


1 

2 


(  9  1 

mr  u/2  - 

.(zmA)1/2  . 


_1 _ 

2  rr 


mr2w/2  -  P. 

V  i- 


x  exp 


(2mr2T)1/2 


1  A  -t 

mm  1 

(2mr2T)1/2 


and  similarly  for  Pg 


Other  Cases : 

Other  cases  may  be  derived  from  these,  if  we 
?2  as  arguments  to  p p g,  p^; . 


P 


Fa  (r,*,Pl,P2) 


etc  . 


r/2) 


add  Pj  , 
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We  then  have, 


Case  D: 


>  P  >  P.  >  mr  w/2  [1 
max  1 


f_8e4J 

1/2  • 

2  2 
Lmr  r  J 

pD  ~  FA  (r  ’  ^  ’ Pmax ’ P2)  +  FB  (r,^’Pl ,Pmax3 


Case  E:  P„  >  P  >  P  >  Pi 

2  max  mm  1 


P  E 


f  =FA  (r ’  *’] "max  ’ P2)  +FB  (r  ’  *>  Pmin  ’  U  +FC  (r  ’  *  >] P1 •  Pmin) 


Case  F:  P  >  Po  >  p  m  >  Pi 

max  2  mm  I 


pF  "  FB  (r  ’  ^  ’  Pmin  '  P2]  +  Fc(r,^,Pl  ,Pmin) 
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